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Abstract 

We consider the evolution of a parton system which is formed in the central region just after a 
relativistic heavy ion collision. The parton consist of mostly gluons, minijets, which are produced 
by elastic scattering between constituent partons of the colliding nuclei. We assume the system can 
be described by a semi-classical Boltzmann transport equation, which we solve by means of the test 
particle Monte-Carlo method including retardation. The partons proliferate via secondary radiative 
gg ggg processes until the thermalization is reached for some assumptions. The extended system 
is thermalized at about t = 1.6 fm/c with T = 570 MeV and stays in equilibrium for about 2 fm/c 
with breaking temperature T = 360 MeV in the rapidity central region. 
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I. INTRODUCTION 



It is generally expected that a Quark-Gluon Plasma(QGP) can be formed after a rela- 
tivistic heavy ion collision even at RHIC energy. However, it has not been clear yet how 
a QGP is formed in a heavy ion collision, as well as what the equation of state is, if it is 
formed. Some of these questions, which are important for the identification of signals of 
the QGP and of the properties of dense strongly interacting matter, are the subjects of our 
study. 

Beginning with the moment of contact of two colliding heavy ions, initially coherent 
quanta, which can be identified as partons, from the projectile nucleus collide with those 
from the target nucleus and can become incoherent particle excitations. They can, as well, 
produce additional particles by radiative processes 0|. There may be several different 
mechanisms responsible for particle production in relativistic heavy ion collision depending 
on the transverse momentum of the produced particles. Those of high transverse momentum, 
which is greater than a threshold momentum po, are usually called minijets and can be 
approximately calculated by the methods of perturbative QCD (pQCD) using the parton 
structure functions of the colliding nuclei. On the other hand, the "soft" partons of low 
transverse momentum, with pt < Po, are thought to be produced by the semi-classical 
decoherence of strong, random gluon fields contained in the colliding nuclei, sometimes 
called color glass condensate 

While the primary partons of high transverse momentum (minijets) are produced near 
the central rapidity region and have relatively high transverse energy, the soft partons are 
distributed over a wide range of rapidities and have lower transverse energy. These soft 
partons can play a role as a thermal bath or background field, which slows down the minijets 
by inducing them to transfer a substantial fraction of their energy to the bath, thus enhancing 
their chance to be thermalized. It was pointed outj^ 7|, however, that the major mechanism 
of energy deposition is through radiative (2 3) processes, suggesting that it is sufficient 
to consider a minijet system in order to explore the chemical and thermal equilibration of 
the system. This is what we are doing here. 

We study the initial phase-space distribution of minijets assuming no correlation be- 
tween momentum and space and another one with correlation, especially those given by 
Cooper at al[lQ], which we will call CMN distribution, in Section II. We briefly describe the 
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semi-classical Boltzmann equations of motion and the Monte-Carlo simulations to solve the 
equations in Section III. We present the results of our numerical solutions and a summary 
in Section IV. 



II. INITIAL DISTRIBUTIONS 

The hard gluons produced after a relativistic heavy ion collision have been studied exten- 
sively by many authors. Its momentum distribution as well as the total number of partons 
can be calculated by the method of pQCD. To see the dominant effect at very early times 
after the onset of a relativistic heavy ion collision, we consider only minijet gluons. Using 
the parton distribution of colliding nucleus A, fi/A{x,Q'^) = fi/N{x,Q^)RA{x,Q'^), where 
fi/N{x,Q'^) is the parton distribution of a free nucleon and Q^) the nuclear ratio 

function, the minijet distribution can be calculated by the 2 — > 2 minijet cross section per 
nucleon in a A + A collision. 



dprdy 



dy2—:^ ^ Xlfi/A{Xl,PT)x2fj/A{x2,PT)^i'j^kl{s, t, u) , (2.1) 
ij, kl 

where K denotes a factor to include the higher order diagrams, which we will set to K = 2 
throughout our study. T{b) is the nuclei geometric factor at an impact parameter b. We use 



111 ] and the EKS98 parametrization 



the GRV98 set of parton distributions for a free nucleon 

n 

for the ratio function Xi and X2 are the Bjorken scaling variables of parton i and j 

in nucleons of two colliding nuclei, respectively, yi and y2 are the rapidities of scattered or 

produced partons. 

The gg scattering cross section at leading order is 

97ra^ til us si 
^™ = ^[3---^--], (2.2) 

where the "hat" symbol on the Mandelstam variables indicates those of a parton. The 
relations between the variables are given by 

xi = PT(e^i +e2'^)/v^, (2.3) 
X2 = prie-y' +e-y')/^/^, (2.4) 
s = X1X2S, (2.5) 

i = -pUi + ey^-y'), (2.6) 

u = -p'^{l + ey'~y-'). (2.7) 



We also rewrite the available kinematic region for convenience 



Po'<PT'<{,r^)\ (2.8) 
2 cosh y 

- log(^ - e-y) <y2< log(^ - e-y), (2.9) 
Pt Pt 

We consider only a head-on collision {b = 0) and use the leading-order expression for the 
running coupling constant = 47r[6o log(Q/AQCD)]~^; where Q is the momentum transfer, 
bo = 11 — |n/, and Aqcd = 200 MeV and nj = 3. The production of minijet gluons strongly 
depends on the transverse momentum cut po. We use po = 1 GeV at RHIC energy, which 
is marginal for the use of pQCD. Note that we set the minimum transverse momentum to 
3e equal to the saturation momentum, where the gluon distribution becomes semi-classical 
^. The total numbers of minijets (integrated over rapidity) is about 4800 at RHIC. These 
numbers are quite large since we set the momentum cutoff for the minijets at the lower limit 
of the range p^ ^ 1 — 2 GeV^ at RHIC. The rapidity distribution is shown in Fig. 1 of 
reference (lOl |. 

To obtain the phase space distribution of minijets, we first assume that there are no corre- 
lation between the momentum and space coordinates of a produced parton. The momentum 
distribution can be given by the minijet distribution, Eq. ()2.H) . 

1 dN^^^ 
p\ dy dpT ' 

We can choose {pT,y) using this distribution and the azimuthal angle with equal weight 
between (0, 27r) by the Monte-Carlo method, by which we select the energy and momentum 
of the produced parton. The spatial distribution of minijets can be deduced from the fact 
that the minijets are produced from elastic scattering between constitutient partons of two 
colliding nuclei. Assuming the space-time point of produced particle is just that of elastic 
scattering, we can estimate the space-time positions based on the classical picture as follows. 
We assume that the probability for having an elastic scattering of constituent partons is 
proportional to the nuclear density overlap of the colliding nuclei. The transverse position 
of elastic scatterings is thus obtained by 

X = r sin 5 cos 0, (2.12) 
y = r sin 6' sin 0, (2.13) 



f{pT,y) = C—-J-J-. (2.11) 



where the distribution of constant nuclear density in the nuclear rest frame is P{r, 6, 0) = 
3rV(47ri?3). 

To obtain the longitudinal position and the time of a collision, we consider a sphere of 
radius R which has a constant density and is moving with a constant velocity v ^ c. We 
assume two identical spheres suffering a head-on collision in the CM frame. At a given 
transverse collision position {x,y), we consider longitudinal tubes of area AA through the 
point {x,y). Because of the Lorentz contraction, the half-length of the tube is D/'j with 



D = B? — — 'xp' and 7 = (1 — f^/c^)^^/^. The tubes from both nuclei begin to overlap 
one another starting at ts = {R — D)/'yv and ending at te = {R + D)/'^v. Note that we set 
the time of first contact between the two colliding nuclei at t = 0. At time t, the overlap 
volume is given by 

AV" = 2AAv{t-ts)h. (2.14) 

The maximun overlap volume is AV = 2AAD/'~f at t = R/'yv. The probability density 
which elastic parton collisions occur is proportional to the overlap volume and the product 
of the two Lorentz contracted parton densities: 

2 2 

m = l^it-ts) (2.15) 

Using this collision time distribution, we select the collision time randomly: 

2£' ^ . ^ 

t = ts + — 2.16 

where ai G [0, 1] is a random number. Once we have t, we can calculate the overlap volume, 
Eq. ()2.14| ). and its longitudinal range, —v(t — tg) < z < +v(t — tg). We again choose the 
longitudinal location z of the collision in this range randomly: 

z = v[2{t-ts)a2-{t-Q], (2.17) 

where 02 G [0, 1] is another random number. 

We also consider the phase-space distributions ()2.H) at RHIC which were recently pro- 
posed by Cooper, Mottola and Navak[lo| (CMN Distribution). The boost noninvariant 
distribution, for example, is given 

nPT,V,y,ro) = C^^^^^^^^j^^,, (2.18) 
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FIG. 1: The transverse momentum distribution of a Monte Carlo sampled parton system as a 
function of transverse momentum in GeV unit. 

where t] is spatial rapidity. is an unknown parameter function, which relates the space- 

rapidity correlation to the transverse momentum and was chosen to be constant (cr^ = 0.28) 

n 

at RHIC energy in Ref. jlO(|. Using this distribution function, we can randomly selects the 
phase-space coordinates (pr, y, v) of ^ach produced parton. We assume in this distribution 
that the parton production time is given by Eq. ()2.16|) . 

Figure shows the transverse momentum distribution of Monte-Carlo sampled partons 
with a minimal transverse momentum po = 1 GeV/c. Most of the produced particles have 
the transverse monenta less than 3 GeV/c and the number of produced gluons sharply 
increases as the transverse momentum gets smaller. 

Figure El shows the rapidity distribution of Monte-Carlo sampled particles, which can be 
compared to the rapidity distribution given in Ref. Q|. The distribution shows that the 
rapidity is flat in the central region, but falls off for \y\ > 3. 

Since we use Monte-Carlo sampling of the distribution, we use several different realizations 
of the initial distribution (events) and average over the events to get the physical values. 
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FIG. 2: The rapidity distribution of the produced partons. 



III. NUMERICAL SOLUTION 



The initial-state gluons, which we will call 



to the Boltzmann equation 



OTimary gluons from now on, evolve according 



p''dMx,p) = I //^W^„[/,(3)/,(4)-/,(l)/,(2)] 

•/ 2 3 4 

/ / / W^.™[/.(3)/.(4)-/,(l)/,(2)] 

J2 Js Ji 



+ 
+ 
+ 




gg^aaa 



imm - /.(i)/.(2)] 



2J3J4J5 




W, 



ifMm) - /.(i)/5(2)] 



2 J3 J A 



1 1 l\[w,,^M?,)f,{^)-w,,^,M^)f,{2)i 

2 3 4 



(3.1) 



where we used the abbreviated notations fg{i) = fg{x,pi;t), q = {u,d, s), q = {u,d,s), and 
J. = J dpi/Ei. The transition rates can be expressed as 



Wgg,^g,g'^ = Sa{S, Q) d^^^P +P,~p'- p[) . 



(3.2) 



Note that we assign the label 1 to the momentum p. We explicitly include the particle 
symmetry factor in the classical limit. The ggg —>■ gg process is not included in the equation 



of motion since the system of high energy partons expands very rapidly so that the available 
volume increases and the channel is suppressed strongly. 

On the other hand, the secondary quarks and antiquarks produced in parton collisions 
satisfy the equations of motion, 



where C = 1/2, if the final state consists of identical particles and C = 1 otherwise. We 
neglect the quantum mechanical effects from Bose enhancement factors (1 + fg) for a gluon 
in the final state and Pauli blocking factors (1 — fg) for final-state quarks or antiquarks. 

We use the parton cascade code, which was developed by us and described in ref. to 
solve the equations of motion for a given initial state. We first initialize all primary partons 
in space-time with their four-momenta and look for the next possible collision. If a collision 
occurs, we select the collision products and their new momenta according to the total and 
differential cross section, respectively, at the space-time of relativistic maximum force point 
between the colliding particles. This procedure is terminated when the evolution time of 
the cascade is greater than the final time set as an input parameter. 

IV. RESULTS AND CONCLUSION 

Fig. El shows the ratio ^ \Px\/Yl \Pz\ ^ small sphere with radius i? = 1.1 fm at the 
center of the coordinate frame as a function of time, where 




(3.3) 




(3.4) 




(4.1) 
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When this ratio becomes unity, an isotropic, equihbrated configuration has been reached. 

As expected, the ratio for the evolution without colhsions remains clearly different from 
unity, which means that there is no kinetic equilibration. The ratio with collisions shows 
the kinetic isotropy from a time of about 8 GeV^^ = 1.6 fm/c to 17 GeV^^ = 3.4 fm/c 
for the system initialized with minijet partons. The simulation includes small angle elastic 
scattering, for which we set the minimum momentum transfer to Pcut = 0.5 GeV/c. We also 
show the results of simulations with hard scattering only, pcut = 1 GeV/c. Even though the 
total number of partons in the system is almost same with 9700, only the system with small 
angle scatterings reaches a kinetically isotropic state. 

On the other hand, the ratio for the system of minijets of CMN distribution always re- 
mains away from unity. We can understand this result as follows: Since the value of the 
spatial rapidity of a produced parton is close to its momentum rapidity, those partons which 
have positive longitudinal momentum, will be produced with positive longitudinal position, 
the right hand side of the collision plane, and the parton with negative longitudinal momen- 
tum will be born at a negative longitudinal position, on the left-hand side of the collision 
plane. This correlation thus makes possible only those collisions between the particles which 
are running in the same longitudinal direction but different transverse momentum. This phe- 
nomena reduces the number of collisions and the collision energy substantially, so that the 
evolution is essentially free streaming. 

Figure lU shows the number density in the sphere as function of time. The density of 
partons for minijets simulations decreases faster than 1/t and from 32 fm~'^ at t = 8 GeV^^ 
to 11 fm^^, which is much higher than the density of valence quarks of a free nucleon, at 
t = 17 GeV~^ in the small sphere at the center. The density equals to the valence quark 
density of a free nucleon at t = 34 GeV~^. Note that we plot the data if the particle density 
of the sphere is greater than that of valence quarks in a nucleon, 1.4 fm~^. 

Figure |31 shows the energy density of the sphere as a function of time, which is almost 
exponential decreasement. The energy density decreases from 54 GeV/fm^ at t = 8 GeV^^ 
to 12 GeV/fm^ at t = 17 GeV^^ for our minijet simulations. 

Figure IHl shows the energy per parton in the sphere as a function of time. The energy 
per parton of free streaming(no collision) is much higher at early stage of expansion than 
those of interacting systems, which means the parton loses a substantial amount of energy 
via scattering. The energy decreases from 1.7 GeV at t = 8 GeV~^ to 1.1 GeV at t = 17 
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FIG. 3: The ratio ^^^"^ in a small test sphere (r = 1.1 fm^) at center region. 
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FIG. 4: The particle density N/V in a small test sphere at central region. 

GeV^^ for a minijet system. If we assume from Fig El that the local system of the sphere at 
the center is thermalized, the relativistic gas has the temperature T = 566 MeV at t = 8 
GeV~^ to 367 MeV at t = 17 GeV~^ for a minijet system and a mixed system, respectively. 
This system is much hotter than the QGP transition temperature ~ 170 MeV. 
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FIG. 5: The energy density E/V of the sphere. 
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FIG. 6: The energy per particle E/N in the sphere. 

Figure m shows the total number of particles in the system. The secondary particles are 
produced within 2 fm/c after the first contact of two colliding nuclei. This shows that the 
chemical freeze-out comes much earlier than the thermal freeze-out even at the center of 
collision. The total number of particles is about 9700 for a minijet system and 5200 for 



11 



o 

CL 



0) 
JO 

E 

o 



10000- 



8000 



6000 



4000 



2000 



10 



Minijets with NO Collision 
Minijets 

Minijets of CMN Distribution 



20 30 

Time(1/GeV) 



40 



50 



FIG. 7: The total number of particles in the system as function of time. 



a minijet system of CMN distribution. As many authors have pointed out, the gg ggg 
process, even though the process has the kinematical constraint, 6{pTXf — coshy)6{y^/2 — 
PTCOshy), which is E^m > 1-54 GeV in our simulations, plays a dominant role in depositing 
energy and thermalization of a system. 

Fig. ISl shows the pseudo rapidity distribution of all of partons in the system at t = 50 
GeV~^. Since partonic collisions are ceased by this time, the distribution will be preserved 
to the final state unless the hadronic interaction changes them. We can see that the ra- 
pidity distribution is narrowed since the high rapidity partons at early time. Fig. El have 
been scattered to produce the low rapidity partons. We also point out that the rapidity 
distribution is different from the final state (hadronic) rapidity distribution in number and 
the width 0- 

We conclude that the system is thermalized even with minijets at RHIC energy if the small 
angle scatterings are included. When we look at the small sphere of radius R = 1.1 fm at 
center, the system is thermalized from t = 1.6 fm/c till about 3.4 fm/c. And the temperature 
is decreasing from 566 to 367 MeV and the particle density from 32 to 11 fm~^. Even after 
this time, the system is still hot untill about 8 fm/c at the center and may cool down further 
after hadronization. The evolution of the system, which has no correlation between position 
and momentum, is very similar in character to the "bottom-up thermalization" scenario of 
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FIG. 8: The pseudo rapidity distribution of parton systems as functions of time. 

Baier et a/.[3|, in which the hard gluons hberated from the colhding nuclei within r ^ Qj^ 
colhde among themselves to produce soft gluons and thermalize via small angle scattering. 
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